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Abstract 

A review of present understanding of the dissipation region in magnetic reconnection 
is presented. The review focuses on results of the thermal inertia-based dissipation 
mechanism but alternative mechanisms are mentioned as well. For the former process, 
a combination of analytical theory and numerical modeling is presented. Furthermore, 
a new relation between the electric field expressions for anti-parallel and guide field 
reconnection is developed. 


1. Introduction 

Magnetic reconnection is often recognized to be the most important plasma transport 
and energy conversion process in space physical plasmas. Magnetic reconnection is a 
likely contributor to the formation and ejection of coronal mass ejecta (e.g., Gosling et 
al., 1995; Antiochos et al., 1999) and to coronal heating (e.g., Forbes and Priest, 1987, 
Cargill and Klimchuk, 1997), and facilitates the entry of solar wind plasma and elec- 
tromagnetic energy into the magnetosphere by either low- or high-latitude magneto- 
pause reconnection (e.g., Paschmann et ah, 1979; Sonnerup et ah, 1981). In the mag- 
netosphere proper magnetic reconnection converts energy stored in the magnetotail 
lobes to plasma internal and kinetic energy. It is also believed to play a role in the 
formation of the auroral acceleration region (Atkinson, 1978; Haerendel, 1987). 
Therefore, magnetic reconnection constitutes a fundamental and ubiquitous element 
of the Sun-Earth connected system. 

Magnetic reconnection relies on the presence of a diffusion region, where collision- 
less or collisional plasma processes facilitate the changes in magnetic connection 
through the generation of dissipative electric fields. This diffusion region is strongly 
localized, extending at most to typical ion Larmor radii. Typical indirect reconnection 
signatures (contrasted with direct observations of the diffusion region) include the 
presence of fast flows and plasma heating associated with magnetic field signatures 
indicative of the establishment of new magnetic connections. Such indirect signatures 
are observed remotely in the solar corona, in the solar wind by the magnetic topology 
of CMEs, and by direct spacecraft observations at the magnetopause and in the mag- 
netotail of the Earth. 


In collisionless plasmas, which we identify with plasmas where collisions are suffi- 
ciently infrequent to provide the necessary dissipation on the time and spatial scales 




consideration, the nature of the dissipation underlying magnetic reconnection 



ter regimes, both in regions of low magnetic field, and in the presence of a magnetic 
guide field. 

In this paper, we present an overview of our present understanding of how particles 
become demagnetized in the diffusion region of the reconnection process. The em- 
phasis of the discussion will be on electrons - however, we point out that many results 
apply equally to ions, after a suitable change of charge and mass. The paper will focus 
primarily on the thermal- or bulk inertia-based dissipation processes that have been 
validated in a large number of numerical models. For completeness, however, we will 
also mention relevant other processes, such as wave-particle interactions that have 
been proposed as candidate dissipation mechanisms as well. 

The structure of the present paper is as follows: Following an introduction of the 
model in the next section, in section 3, we will present a revised theory of guide-field 
magnetic reconnection. Section 4 presents an analysis of our present knowledge of 
anti-parallel magnetic reconnection. Both sections will include su mm aries of transla- 
tionally invariant and fully three-dimensional models. Finally, section 5 will feature a 
summary and a presentation of open questions that pertain to diffusion region physics. 


2. The model and modeled configuration 

Most of the results presented in this overview are based on the application of particle- 
in-cell codes. The system studied by many researchers is similar to that selected in the 
GEM reconnection challenge (Bim et al., 2001). This system is described in the fol- 
lowing. 

During the discussions, we will employ dimensionless quantities. For this purpose, we 
normalize densities by a typical density no in the current sheet, the magnetic field by 
the asymptotic value of the reconnection magnetic field Bo . Ions are assumed to be 
protons (mass m p ) throughout, and length scales by the ion inertial length c/co,, 

where the ion plasma frequency co , = ^e 2 n f) / £ 0 m p is evaluated for the reference den- 
sity. Velocities are measured in units of the ion Alfven velocity 
v A = B 0 / ^j/u (] m p n () based on the reference magnitudes of magnetic field and density. 
The electric field is measured in units of E 0 =v A B 0 , the pressure in units of 
p 0 = B 0 ! jj . 0 , and the current density is normalized to j 0 = (O i B 0 / cp 0 . 


The poloidal magnetic field, a modified Harris sheet (Harris, 1962) with a current 
sheet half width of half the ion inertial length, is of the following form: 

B x = tanh(2z)+a 0 ^ / Z__ cos ( Ittx/LJ sin ( nz/L. ) (4a) 

B : = -a 0 2 jt/L z sinf 2 m/L x ) cos ( m/L, ) (4b) 

The perturbation amplitude ao=0.1 leads to an initial value of the normal magnetic 
field of about 3% of Bo. Similar to earlier studies (Hesse et al., 2002), we here employ 
a constant magnetic field component directed along the main current flow 

By = B y o (5) 


where B y0 is the initial value of the guide field magnitude. The choice of ion-electron 
mass ratio varied by study; it ranges from m t !m e =25. in the original GEM chal- 


lenge problem to m, Jm e =1836 in implicit particle-in-cell calculations (Ricci et al.. 









focus on one representative timestep, t=16. The magnetic field and ^/-directed current 
density are shown in Figure 1 . The initial, X-type perturbation leads to a reconfigura- 
tion. The most prominent difference to anti-parallel merging is the inclination of the 
reconnecting current sheet with respect to the X-axis. The figure also indicates the 
presence of a very thin current sheet in the central reconnection region, which is likely 
associated with electron demagnetization. 

The most interesting locations for the calculation of the electron distributions are 
those where the electron pressure anisotropy provides the reconnection electric field. 
Figure 2 shows the relevant nongyrotropic components of the tensor, obtained from 

( 8 ) 

with 

+ (») 

where enq. (9) represents the gyrotropic component to the full electron pressure ten- 
sor. Figure 2 shows a strong localization of the pressure nongyrotropy near the current 
sheet and magnetic field reversal in the green region of fig. 1 . 

The thinness of the current sheet is demonstrated by Figure 3, which displays a cut in 
the z direction across the diffusion region current sheet. The figure shows the electron 
current density, and the nongyrotropic contribution to the pressure tensor component 
P y:e . A comparison with the magnitude of the local collisionless skin depth, shown at 
the bottom, indicates that the spatial scales of both current and pressure anisotropy are 
smaller. We will show below that the relevant scale here is indeed given by the local 
electron Larmor radius. 

In order to shed light on the underlying physical process, we now investigate particle 
distribution functions. For this purpose, we collect particle distributions at the location 
13.1 <*<13. 2 and -0.0 l<z<0. 04 just above the reconnection region, where the nongy- 
rotropic y-z pressure tensor component provides the gradient, which supports most of 
the reconnection electric field. 

The relevant reduced distribution F(v y ,v z ) is displayed in Figure 4. Here, the green ar- 
row indicates the magnetic field direction, which is essentially exclusively in the y 
direction (B : is negligible). In addition, the white ellipse is drawn to guide the eye in 
the explanation of the nongyrotropy. Specifically, we point out that there is slightly 
more phase space density for positive v_- than for negative v : if v y is large and negative, 
and the opposite holds for positive v y . This slight asymmetry indicates a preference 
for more acceleration in the negative v y direction if particles move away from the cen- 
tral current sheet, and it is sufficient to explain the nongyrotropy of P y:e . 

Preferential acceleration for particles closer to the center of the current sheet is due to 
the strong localization of the parallel electric field, shown in Figure 5. The region of 
particle collection is marked by the black rectangle. The (electron Larmor) scale 
length of the electric field variation implies that particles with instantaneous gyrocen- 
ters below the observation point will have been accelerated more than those with in- 
stantaneous gyrocenters above the observation point. Accordingly, on average, parti- 
cles with v-<0 will have been less accelerated in the negative y direction than particles 
with v->0. This idea is captured in Figure 6* which, show# a disteifeution inaction 













pressure tensor becomes: 


> 1 dQxyze 

yze n ey dz 


( 13 ) 


Further progress therefore requires an analytic expression of an evolution equation for 
the entire heat flux tensor. The heat flux tensor is defined in the electron center-of- 
mass system as: 


Q = m s J d 3 n(u - v){u - v)(u - v)f s 


(14) 


Here, f s and m s denote the distribution function and mass of plasma species s, u the 
phase space velocity, and v the bulk flow speed. An evolution equation for Q is ob- 
tained by multiplying the Vlasov equation by uuu and integrating over phase space. 
The result needs to be transformed into the center-of-mass frame of species s in order 
to derive an evolution equation for Q. After a considerable amount of algebra, one 
finds for the components of the electron heat flux tensor (index e omitted for simplici- 
ty) (e.g., Goswami et ah, 2005) 
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Here % is the usual, totally antisymmetric tensor. We note that the expression quoted 
in (Hesse et al., 2004) was incorrect. We here redevelop the subsequent theory below 
and we will show that the resulting approximate expression for the reconnection elec- 
tric field (Hesse, 2006) remains unchanged. 


Equation (15) relates the time evolution of Qjk to lower order moments such as pres- 
sure and velocities, as well as to the fourth order tensor T ijkl . The last term in (15) rep- 
resents the effects of particle cyclotron motion on the heat flux tensor. Clearly, (15) is 
invariant under change of order of indices, leading to a totally symmetric heat flux 
tensor. 


Further progress toward a simple scaling relation requires simplifying assumptions. 
Neglecting time dependence and the 4-tensor, an expression for Q xyz can be obtained 
from the xv\x-component of (15) 



( 25 ) 


The next step is to look at P xy . The analogue to eqn. (13) is: 
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The corresponding electric field becomes, in leading order: 
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The second term vanishes by symmetry. We therefore get: 
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Using a similar argument as above the tensor component becomes again: 
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The first derivative becomes, in good approximation: 
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The second derivative becomes: 

l 2 

dx 


i/ 

"> 2 *2xyz 


P dv d 2 v v 

xr x y_ 

Q, dx dx 2 


(28) 


(29) 


(30) 


where we have used the symmetries to eliminate some derivatives. The electric field 
contribution becomes, from (27): 

P dv d 2 v 

“ x " (31) 
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which leads to the final result: 
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Taking (28) and (35) together and assuming a weak density variation (consistent with 
the incompressibility assumption) yields: 

„ 1 t dv.. 


2 en L dx 

which is the same result as Hesse (2006). 
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(33) 


Eqn. (33) offers interesting insight into the physics of the underlying reconnection 
process. Based on distribution function analysis we argued above that accelerated par- 
ticles preferentially leave the reconnection region, whereas particles, which have, on 
average, been accelerated less tend to enter. In the absence of acceleration of the re- 
connection electric field the current density would therefore decay - which validates 
our gedankenexperiment above. 


Comparison with (12) reveals that the coefficient K (analogous to a diffusion coeffi- 
cient) is: 


K = 


L r 2<K 

2 L dx 


(34) 


This expression shows the following. First of all, the square of the Larmor radius in- 
dicates that the electric field scales with the square of the ratio of critical length scale 

and Larmor radius in the guide field. Clearly, this implies that the actual current layer 




tion was first described by Horiuchi and Sato (1994). Hesse and Winske (1998) per- 
formed a particle-in-cell simulation of collisionless reconnection in the GEM geome- 
try'. Their analysis showed indeed the presence of electron pressure nongyrotropies 
near the X point. The magnitude of these tensor components proved to be sufficient to 
provide the reconnection electric field via the expression (3). 

These studies were taken a step further by Hesse et al. (1999), who investigated the 
effect of different electron masses on the collisionless dissipation process in the re- 
connection region. The target of the investigation was to study whether different 
physics in the diffusion region might lead to different dissipation, thereby influencing 
and potentially changing the larger scale behavior of the system under investigation. 

In order to test the dependence of the reconnection process on the assumed electron 
mass, Hesse et al. (1999) performed a set of simulations with varying electron masses, 
ranging from mj m e =9 to m, l m e = 100. 

We here discuss a more recent calculation, which employs m i /m e = 100. The system 
size is Z^=100 and L z =50, and a total of 10 9 particles. Figure 8 shows the magnetic 
field evolution with the current density color-coded. Arrows indicate the locations, 
where particles are collected to form distribution function (see below). The figure 
shows that magnetic reconnection, initiated by the initial perturbation, causes large 
changes to the magnetic field and current density distribution. The figure demon- 
strates that the current density exhibits a saddle-point at the location of the reconnec- 
tion region. The latter feature becomes most prominent at later times. 

The relevant off-diagonal components of the pressure tensor are shown in Figure 9. 
Similar results were found also by many other investigations, e.g., Pritchett (2001a). It 
is apparent from the discussion above that the processes responsible for these electron 
pressures rely on the inertia of individual electrons, which contributes to all of the flu- 
id terms on the RHS of eqn. (2). 

Figure 10 shows reduced distribution functions, similar to those presented above for 
finite guide field, at the location indicated by arrows in fig. 8. The left panel shows 
F(v z ,v y ) in the inflow region, and F(v x ,v y ) is shown in the outflow region. Qualitative- 
ly, the figure shows the same features as the ones found in the guide field model 
above: On average, more hot and accelerated plasma is leaving the diffusion region 
than is entering. We will come back to this point below. 

In order to estimate the magnitude of the reconnection electric field, Hesse et al. 
(1999) realized that the scale of the electron diffusion region is given by electron orbit 
excursion in a field reversal, the so-called “bounce width” (6). A corresponding scale 
can be derived for the bounce motion between the downstream magnetic field compo- 
nents B z . The approximate bounce orbit is sketched in Fig. 11. This result can be used 
for an estimate of the electric field from (3): 
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Here the values of the pressure tensors are to be taken at the edges of the current 
sheet, where electrons begin to become magnetized. In these regions, the pressure ten- 
sor components can be approximated,by,(Kuznetso\a et al., 1998): 
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We approximate the momentum density derivative in the following way: 
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where we have neglected the electron flow velocity at the edge of the current layer 
and assumed that the dominant contribution to the current is by the z derivative of B x . 
Assuming L~b as defined in (6) we get: 
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Comparing (42) and the RHS of (40) leads to: 


ry m f «A,v,max 


\m n v 

\ e e ey 


„ 1 rr — 0V 1 1 , 3v , 

It 




the desired expression for While we emphasize the approximate nature of the 
present derivation, this analysis shows that the electric field expressions (33) and (38) 
are of essentially the same nature. 


Thus, the diffusion region in anti-parallel reconnection appears fairly well understood. 
However, there are other modes and instabilities that may change the reconnecting 
system substantially, without necessarily changing the bounce motion-based dissipa- 
tion process. For example, the kinetic drift-kink instability (DKI) (Zhu and Winglee, 
1996, Pritchett et al., 1996, Ozaki et ah, 1996) has been suggested to substantially 
modify the structure of the reconnecting current sheet. Although analytical theory 
(Daughton, 1998, 1999) and numerical models (Hesse and Bim, 2000; Zeiler et ah, 
2002) show that DKI growth rates appear to be small for realistic ion-electron mass 
ratios, there remains some debate whether the DKI or the related sausage mode 
(Buechner and Kuska, 1997; Pritchett et al, 1996; Silin and Buechner, 2003) may play 
a role in furthering the onset of magnetic reconnection. 


Finally, laboratory experiments have made great progress toward the objective to pro- 
duce truly collisionless magnetic reconnection in laboratory devices. Recent studies 
have shown that, although (39) is reproduced in particle-in-cell simulations (Dorfman 
et al., 2008), laboratory measurements fail to reproduce a sufficiently strong electric 
field based on the same equation. This feature may be due to effects not present in the 
simulations models, or they may be an artifact of residual collisionality in the experi- 
ment. 


5. Summary 

The ability to describe the diffusion region is fundamental to our understanding of 
magnetic reconnection. Consequently, the dissipation region physics has attracted the 
attention of a substantial number of investigations. Early studies recognized that, on a 
microscopic level, electron orbits in the diffusion region should differ substantially 
from gyro-drift motions. Theoretical studies of Ohm’s law in the diffusion region in- 
dicated that collisionless processes that provide reconnection electric fields should be 









which, when sufficiently strong, contributes substantially to the overall pressure bal- 
ance. For stronger guide fields, less plasma heating is therefore required to maintain 
pressure balance. This effect is evident when comparing figs. (4), (7), and (10). As 
expected, we find comparatively less (perpendicular) heating in the guide field model 
than in the anti-parallel case, where there is no guide field available to provide mag- 
netic pressure in the diffusion region. The properties above also apply in asymmetric 
magnetic reconnection, even though the overall morphology may be more complex 
than in the symmetric cases discussed here (Pritchett and Mozer, 2009). 

In summary, more research is required to provide closure to the diffusion region phys- 
ics problem. It appears essential that this research be conducted by means of tailored 
space investigations such as those offered by NASA’s Magnetospheric MultiScale 
(MMS) mission, and that this research be challenged by close interactions with labor- 
atory experiments. While the ultimate answer remains elusive today, it is clear that 
this pursuit will be one of the most exciting endeavors in modem plasma physics. 
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Figure 4. Reduced distribution F(v y ,v z ). The distribution exhibits a small but signifi- 
cant nongyrotropy about the magnetic field direction green arrow, which is revealed 
by comparing to the white ellipse. The arrows indicate regions of enhanced electron 
flux. 

B, ve and E_par, t= 16.02 * 10 
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Figure 5. Parallel electric field (color), magnetic field lines (white), poloidal electron 
bulk flow (vectors), and particle collection box. 



Figure 6. Illustration of the generation mechanism on the asymmetry in F(v y ,v : ). Parti- 
cles with gyrocenters below (in z) the dissipation region have experienced, on aver- 
age, more acceleration in the negative y direction than particles with higher gyrocen- 
ters. The result is a nongyrotropy of the type shown in Fig. 4. 










Reconnection electric field : Acceleration, heating, flux transport 


Figure 12. Schematic of the role of the reconnection electric field. 
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